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We investigate the fluctuations of the condensate in the 
ideal and weakly interacting Bose gases confined in a box of 
volume V within canonical ensemble. Canonical ensemble is 
developed to describe the behavior of the fluctuations when 
different methods of approximation to the weakly interacting 
Bose gases are used. Research shows that the fluctuations of 
the condensate exhibit anomalous behavior for the interacting 
Bose gas confined in a box. 
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I. INTRODUCTION 

The experimental achievement of Bose-Einstein con- 
densation (BEG) in dilute alkali atoms spin-polarized 
hydrogen Q and recently in metastable helium has 
enormously stimulated the theoretical research 0] on the 
ultracold bosons. In particular, fluctuations (S^Nq) of 
the mean ground state occupation number Nq have been 
recently thoroughly investigated in a series of papers. 
Apart from the intrinsic theoretical interest, it is foresee- 
able that such fluctuations will become experimentally 
testable in the near future §. 

It is well known that within grand canonical ensemble 
the fluctuations of the condensate are given by <^(5^iVo) = 
-^0 (-^0 + 1) ^ V'^, implying that 6Nq becomes of order 
N when the temperature approaches zero. To avoid this 
sort of unphysically large condensate fluctuations, canon- 
ical (or microcanonical) ensemble has to be used to inves- 
tigate the fluctuations of the condensate. Within micro- 
canonical and canonical ensemble, the fluctuations of the 
condensate have been studied in a systematic way in the 
case of the ideal Bose gas |6|-|l3| . Recently, the question 
of how interatomic interactions affect the fluctuations of 
the condensate has been the object of several theoreti- 
cal investigations [^-19|. Giorgini et al. found the 
anomalous behavior of the fluctuations in a weakly inter- 
acting Bose gas confined in a box within the traditional 
particle-number-nonconserving Bogoliubov approach. In 

Pthe fluctuations of the condensate follows the law 
A^o) ^ V^l'^ . However, Idziaszek et al. considered 
that the fluctuations are proportional to the volume. Re- 
cently, Kocharovsky et al. |l8[ supported and extended 



Although the correction to the ground state occupa- 
tion number due to interatomic interaction has been 
clearly discussed within grand canonical ensemble |2^] 
and canonical ensemble , the role of interaction on the 
condensate fluctuations of the weakly interacting Bose 
gas is still an open and unsolved problem. Different from 
the ground state occupation number, different models of 
describing the weakly interacting Bose gas will lead to 
vastly different prediction concerning the fluctuations of 
the condensate. 

The purpose of this paper is to present a unifled 
method to calculate the fluctuations of the condensate 
when different ways of approximation to the weakly in- 
teracting Bose gases are used. Within the canonical en- 
semble we give the distribution function of the ground 
state occupation number for the ideal and interacting 
Boson system in a box. We obtain the fluctuations of 
the condensate from the distribution function. In par- 
ticular, we found that the distribution function is not 
Gaussian function in the case of interacting Boson sys- 
tem in a box. The paper is organized as follows. Sec. II is 
devoted to outline the canonical ensemble, which is de- 
veloped to discuss the fluctuations of the condensate for 
the ideal Bose gas in a box. In Sec. Ill we investigate the 
fluctuations of the interacting Bose gas based on the low- 
est order perturbation theory. In Sec. IV the fluctuations 
are calculated based on the Bogoliubov theory. Finally, 
we give a discussion and summary of the results in Sec.V. 



II. MEAN GROUND STATE OCCUPATION 
NUMBER AND FLUCTUATIONS IN THE IDEAL 
BOSE GASES 

Let us start our investigation on the fluctuations of the 
ideal Bose gases in the frame of canonical ensemble. Ac- 
cording to the canonical ensemble the partition function 
of the N non-interacting bosons in a box is given by 
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exp [-/3{S7V„£„}], 



(1) 



where Nn and £„ are the occupation numbers and energy 
level of the state n = {nx,ny,nz} respectively. (3 = 
l/ksT. In (|l|) the energy level of the system takes the 
form 



the results of the work of Giorgini et al. |14| using the 
particle-number-conserving operator formalism. 



n^n^ (r 



2mL^ 



(2) 



1 



Separating out the ground state n = from the state 
n 7^ 0, we have 



N 



[N] = J2 {e^P i-P^oeo] Zo [N - No)} , (3) 

No=0 

where Zq (N — Nq) stands for the partition function of 
a fictitious system comprising N — Nq non-interacting 
bosons. Assuming Aq {N — Nq) stands for the free en- 
ergy of the fictitious system, 

Ao{N-No) = -ksT In Zq (N - No) (4) 

From (||) and the partition function Zideai [N] be- 
comes 
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Z^eal [N] = J2 ^O) 



(5) 
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where q {N, No) = -/JA^oEo - /3Ao [N - No). Obviously 
exp [q{N, No)] /Zideai [N] represents the probabihty to 
find No atoms in the condensate. We will give the distri- 
bution function of the ground state occupation number 
in the following. 

Let us first investigate the largest term in the sum 
of the partition function Zideai [N]. Assume the num- 
ber of the condensed atoms is Nq in the largest term 
of the partition function Zideai [N]. The largest term 
Zq {N — Nq) is determined by the requirement that 
glr^q{N,No)\No=NS=0. «e. 



/3eo - P^Ao (N -NP)=0 



(6) 



The calculations of the free energy (A^ ~ A^o ) non- 
trivial because there is a requirement that the number of 
the particles is A^ — A^q in the summation of the partition 
function Zq {N — Nq). Using the saddle-point method 
developed by Darwin and Fowler ||2^ it is straightfor- 
ward to obtain the free energy (A^ ~ A^o ) of the ficti- 
tious system. 

Ao {N ~ N^) = {N- N^) ksTlnz^Q V^g,,^ (zg) , 

(7) 



where A = ■^2iT(ih^/ m is the thermal wavelength. Zq is 
the fugacity of the A^ — A^q non-interacting bosons and is 
determined by the equation 



N-N^o=T. 
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f;^ exp [e„/fcsr] (^S) 



= TI53/2 (^o) • 



(8) 



From and (||) one finds 



P^^Ao{N^N^o)=^^4- 



(9) 



Combining and (^) one obtains In Zq — fiso . There- 
fore the most probable value A^q is determined by 



N^o^N-Y, 



exp [(e„ - eo) /ksT] - 1' 



(10) 



A"q is exactly the mean occupation number of the conden- 
sate atoms in the frame of the grand canonical ensemble. 
For sufhciently large A^, the sum X]jVo=o '^"^ (H) maybe 
replaced by the largest term, for the error omitted in do- 
ing so will be statistically negligible. In this case, ( |l0| ) 
shows the equivalence between canonical ensemble and 
grand canonical ensemble for large A^. From (|l^), below 
the critical temperature Nq is determined by 



Nl = N-^(:(im = NU 



3/2N 



(11) 



where T\ 



•0 
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m/2)]V^ ^ ivT^^ is the transition tem- 
perature of the ideal Bose gas. 

Other terms in (||) will contribute to the fluctuations 
of the condensate, and lead to the deviation of the mean 
occupation number (A^o) from the most probable value 
NP. When iVo ^ gf^g (A^, A^o) 0. Assuming 
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qiN,No)=a{N,No), 



(12) 



and repeating the saddle-point method, one obtains 

1 



Nq = n-J2 



exp [(£„ - eo) /keT] exp [-a (A^, A^o)] - 1 ' 



From (|ll|) and (13) one gets 



a{N,No) = - 



X' iNo-NS)\ 



1/2 
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where 6 {No - A^^) is a Sign function. {No ~ N^) = 1 
when A^o > A^^, and {No ~ N^) = -1 when A^o < 
NP. To obtain (|l|) we have used the expansions 
53/2 (1 ~ <5) ~ C(3/2) — 2\/ttS and the approxima- 
tion exp[-a{N,No)] w 1 - Q!(A^,A^o). From (|l2|) and 
( p^ ) one obtains easily the following result for q {N, No). 



i-No 

q {N, No)= a {N, No) dNo + q {N, N^) 
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\No-NPQt + q{N,NPQ). 
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The partition function Zideai [N] is thus 



(15) 
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The coefficient in (||) is A = 21 -k^ x En^ol/n" = 0.105. 
The dashed Une is larger than our result because of the 
approximation in iQ. In |14| (j^A^o) = X^ht^o /n: where 
/n — (exp {sn/kBT) — 1) ^. For the convenience of cal- 
culations, /n is approximated as en/ksT for low energy 
atoms. However, this approximation is used also for the 
atoms whose energy level is larger than ksT. Obviously, 
this approximation will lead to the fluctuations become 
larger. On the other hand, (|2^) holds in the canonical 
ensemble except near and above Tc, while our analysis 
holds also for the temperature near T°. In Fig. 2 the dot- 
ted line shows the numerical result of Wilkens et al. . 

In Fig. 2 the arrow marks the temperature which 
corresponds to the maximum fluctuations (S'^Nq) 
Below the temperature T^,, from (|19|), one obtains the 
analytical result for the fluctuations of the condensate. 



(22) 



It is interested to find that the coefficient differs by a 
factor 2, compared to (|2|). The thin solid line shows 
(H) in Fig.2. 

We should note that our results are reliable although 
the disputable saddle-point method is used to investigate 
the fluctuations of the condensate. It is well known that 
the applicability of the saddle-point approximation for 
the condensed Bose gases has been the subject of a long 
debate Recently, the analysis in showed that 

the fluctuations are overestimated, and do not appear to 
vanish properly with temperature using the usual saddle- 
point method. Our discussions for fluctuations are rea- 
sonable because of two reasons, (i) As proved in 
the free energy (Q) of the non-interacting Bose gases is 
still correct, even when carefully deal with the failure of 
the standard saddle-point method below the critical tem- 
perature, (ii) In the usual statistical method (A'o) and 
(j^A'o) are obtained through the first and second partial 
derivative of the partition function respectively. When 
saddle-point approximation is used to calculate the par- 
tition function of the system, the error will be overesti- 
mated in the second partial derivative of the partition 
function so that we can not obtain correct fluctuations 
of the condensate in the usual method. However, in this 
paper what we used is the reliable result (^. The distri- 
bution function of the ground state occupation number 
is obtained directly from (^, without resorting to the 
second partial derivative of the partition function. (A'o) 
and {6'^ No) are obtained from the distribution function 
in this paper. 



III. FLUCTUATIONS OF THE CONDENSATE 
BASED ON THE LOWEST ORDER 
PERTURBATION THEORY 

In the case of interacting Bose gases, the role of inter- 
actions on the fluctuations of the condensate is still an 
open and unsolved problem. Giorgini et al. p^ ] predicted 
the anomalous behavior of the fluctuations in a weakly 
interacting Bose gas confined in a box, while Idziaszek 
et al. considered that the fluctuations are normal. 
Researches show that different model of approximation 
to the interacting Bose gases will lead to different pre- 
dictions concerning the fluctuations of the condensate. 
The method developed to obtain the fluctuations of the 
ideal Bose gas in this paper can be used straightforwardly 
to discuss the fluctuations of the interacting Bose gases 
when different models of approximation are adopted. 

Let us first discuss the fluctuations of the condensate in 
the case of the lowest order perturbation theory, which 
is also discussed in In terms of the lowest order 

perturbation theory, the partition function of the system 
within the canonical ensemble is given by 



[N] 



exp[-/3(S7V„e„-f £;„t)], (23) 



where the interaction energy of the system takes the form 
p5l,p6| 



Eint — 



mV 



(24) 



In (^^ a is the scattering length. Separating out the 
ground state n = from the state n 7^ 0, one obtains 
the following form for the partition function 



N 



Z^nt [N] = J2 {^^P hf^^oeo - f3E,nt] Zq [N - A^o)} , 

(25) 



Nn=0 



where Zq (N — No) stands for the partition function of a 
fictitious A^ — A'o non- interacting bosons. Using the free 
energy Aq {N — A'o) of the fictitious system, the partition 
function is thus 



N 



Zint 

[N] = [q {N, No) 



(26) 



No=0 



where q {N, No) takes the form 

q [N, No) = -l3Noeo - PE„u - pAo {N - No) . (27) 

Analogously to the case of ideal Bose gases, let us 
first investigate the largest term in the sum of Z^t [N]. 
The largest term is determined by the requirement 
-^f^qiN,No)\isig=N'' = 0- Therefore one obtains the 
most probable value A^q of the interacting Bosons. 
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easy to recover the fluctuations of the ideal Bose gases. 
According to (39) the leading contributions to the fluc- 
tuations are anomalous, ie. proportional to N^/^, while 
there are also normal contributions proportional to N 
due to interatomic interaction. This conclusion contra- 
dicts with that of which predicts normal behavior, 
where the lowest order perturbation theory is also used 
to discuss the fluctuations of the condensate. 



is the occupation number of the quasi-particles. From 
the form of /„, we can construct the partition function 
of the quasi-particles in the frame of canonical ensemble. 



From 



Zb = exp 

{n} 

Zb becomes 



(46) 



IV. FLUCTUATIONS OF THE CONDENSATE 
BASED ON BOGOLIUBOV THEORY 

Let us investigate the fluctuations of the condensate in 
the framework of Bogoliubov theory of a uniform weakly 
interacting Bose gas confined in a box. According to 
Bogoliubov theory |^,^ , the total number of particles 
out of the condensate is given by 



(40) 



where 



2 I 2 



9 N 1/2 



2iF 



(41) 



Zb^ E 

{T.Nn = N} 



-/^E^nf 
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where e*^^ — ^nl (^n + ''^n) can be regarded as an effec- 
tive energy level of the thermal atoms. In this case Zb is 
the partition function of a fictitious Boson system com- 
prising N non-interacting Bosons whose energy level is 
determined by e'^^ . From (|4^) the most probable value 
Nq is given by 



n#o exp 



/kBT 



1 



(48) 



Obviously the occupation number of low |n| in ( p8| ) co- 
incides with that of (BO). Analogously, other A^'o is thus 



grip 



(42) 



and /n is the number of quasi-particles present in the 
system at the thermal equilibrium. 



1 



exp [e^/kBT] ~ 1 



(43) 



In addition, the energy of the quasi-particles entering 
(|l|) and (|4|) is given by the well known Bogoliubov spec- 
trum 



;f = ((en + gno f - g^nfi^ ' , 



(44) 



where g — AttTi a/m is the coupling constant, and 
no = Nq/V is the condensate density. At low |n| = 

^Jn"^ + ny + n1, one obtains ~ cx l/|n| and 

/n l/|n|. This results in l/|n|^ divergence in (^) 
at low |n|. Although this sort of divergence will not lead 
to large contribution to the number of low energy quasi- 
particles, it gives leading contribution to the fluctuations 
of the condensate, as pointed out in [Q. We will inves- 
tigate the fluctuations due to low energy quasi-particles 
in the following. 

In (^ ) Nn can be regarded as the effective occupation 
number of the thermal atoms, while 



= /n 



(45) 



A^o = A^- 



E 



n^o exp 



[et," ~ el"^ /^bT] exp [-a {N, No)] - 1 ' 



From ( [48| ) and (49) one gets 

A^o - A^o^ 



a{N,No) 



(49) 



(50) 



where we have used the approximation /„ sa kBT/e^ 
for low energy quasi-particles. Therefore the Gaussian 
distribution function of the system is given by 



Gb {N, No) = exp 



(No-K)' 



exp 



2En^0 i< + vl) fll 

(C(3/2))"/^ (A^o-A^g) 
(27r)^ AA4/3t2 



(51) 



Obviously the mean occupation number (Aq) and fluc- 
tuations (J^A'o) is given by 



Eno=oNoGb (N^Nq) 
Eno^oGb {N,No) 



{No) 



(52) 
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